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Abstract 

The effect of partial absorption on a diffusive particle which stochas- 
tically resets its position with a finite rate r is considered. The particle 
is absorbed by a target at the origin with absorption 'velocity' a; as 
the velocity a approaches oo the absorption property of the target 
approaches that of a perfectly- absorbing target. The effect of partial 
absorption on first-passage time problems is studied, in particular, it 
is shown that the mean time to absorption (MTA) is increased by an 
additive term proportional to 1/a. The results are extended to multi- 
particle systems where independent searchers, initially uniformly dis- 
tributed with a given density, look for a single immobile target. It is 
found that the average survival probability P'^'" is modified by a mul- 
tiplicative factor which is a function of 1/a, whereas the decay rate of 
the typical survival probability P^^^ is decreased by an additive term 
proportional to 1/a. 

PACS numbers: 05.40.-a, 02.50.-r, 87.23. Ge 

1 Introduction 

Strategies to tackle search problems typically involve a mixture of local steps 
and long-range moves[T]. The local part of the search is often considered as 
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a diffusive process whereas the long-range moves may be drawn from some 
Levy distribution [2] or, as studied recently, may be the stochastic process 
of 'resetting' the search to some preferred position[21 SI El E]- From a the- 
oretical perspective the dynamics of the searcher seeking its target may be 
simply modelled by a diffusive particle or random walker which is absorbed 
by the target when it arrives at the target's position. Then a measure of the 
efficiency of the search is the mean time to absorption of the process [7]. 

One simplification inherent in this description is that absorption occurs 
instantaneously i.e. the searcher instantaneously finds or identifies the target 
upon contact. Of course real searches are less reliable and other factors come 
into play e.g. human or animal error may occur in foraging or searching 
for a face in a crowd; there may be stochasticity in the biochemical binding 
reaction of a protein searching for a target promoter site on a length of DNA. 
Moreover, in certain situations the target itself may switch stochastically 
between states that are visible or invisible to the searcher, for example a 
fleeing prey hiding from a predator or a binding site that switches from 
being available to unavailable through external factors. 

In the reaction kinetics literature a partial reaction is often modelled by 
a so-called ^radiation boundary condition^ where the probability density flux 
into the target site is proportional to the probability density at the target 
site[7l El E]. However, here we find it more convenient to instead use a 'sink' 
term in the master equation which represents the loss of probability density 
from the target site. It has been shown flUi [TT] that the two approaches are 
equivalent. 

Our aim in this paper is to investigate the effects of a stochastic absorption 
process of the searcher by the target. We shall consider the diffusion process 
with stochastic resetting and introduce a finite rather than infinite absorption 
rate of the searcher by the target. We will refer to this as partial absorption 
as opposed to full absorption. The mathematical problem then becomes, 
rather than computing the mean first passage time to a target site, the mean 
time to absorption by that site. We investigate the robustness of the results 
of [21 H] to partial absorption. We find that generally for high absorption 
rate the changes are small, so that the results are indeed robust. However, 
the changes in some properties are more significant than others. For example 
the mean time to absorption differs from the mean first passage time by an 
additive constant whereas the survival probability of the target decays with 
time exponentially with a rate that depends on the absorption rate. 

The paper is organised as follows. In section 2 we define the model, 
discuss boundary conditions relevant to partial absorption and identify im- 
portant dimensionless quantities in the problem. In section 3 we compute 
the survival probability of a single diffusive particle and compute the asymp- 
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totics. In section 4 we study the mean time to absorption and how this may 
be minimised with respect to the resetting rate. In section 5 we consider 
many independent particles and compute the average and typical survival 
probabilities for a target at the origin. We conclude in section 6. 



2 Model 

The model considered is that of a diffusive particle which stochastically resets 
its position and can be absorbed at the origin. To study the probability 
distribution, p{x,t\xQ), of the particle at time t after having started from 
initial position xq we construct the (forward) Master equation 

dp{x,t\xo) d'^p{x,t\xo) / .1 N , A (n +\Jif \ fT\ 
— = D — rp{x,t\xo} +r6{x-xo} - ap{0,t}6{x} . (1) 

The final term in the Master equation represents the partial absorption pro- 
cess at the origin. The particle is absorbed with rate a6{z), where the con- 
stant a has dimensions of velocity and we refer to it as the absorption velocity. 
The first term on the right hand side corresponds to diffusion with diffusion 
constant D, and the second and third terms correspond to resetting from all 
points X to the point Xq with rate r. 

We study the survival probability q{z,t) of a diffusive searcher (particle) 
after time t given that it started at position z. It is convenient to consider 
the backward Master equation for the survival probability which reads as 
follows [3]: 

= D^^^ - rqiz, t) + rg(xo, t) - aqiO, t)5(z) . (2) 

The process of resetting to position xq corresponds to the two terms on the 
right hand side proportional to the resetting rate r, which represent loss from 
all positions and gain at Xq of probability respectively. The final term on the 
right hand side represents absorption at the origin. The system with no 
absorbing target site (or a target which does not interact with the searcher) 
is recovered by taking a = 0. The total absorption limit can be recovered by 
taking a — )■ oo, as we now discuss. 



2.1 Boundary Conditions 

The most commonly studied boundary condition in the diffusive process lit- 
erature is the totally absorbing boundary condition wherein, once the particle 
meets the boundary, it is certain that the particle is absorbed and leaves the 
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system instantaneously. The totally absorbing boundary condition is there- 
fore defined, in terms of the probability density p{x, t) of the particle at the 
boundary xb, as p{xB,t) = 0. 

To describe a system which includes a 'refiecting wall', it is natural to use 
what is commonly referred to as the ^reflecting boundary condition'. Here, 
because all particles which try to reach the boundary location are 'refiected' 
back away from the boundary, the net fiux into the boundary site is zero. 
This means that this boundary condition can be straightforwardly expressed 
as 

dp{x, t) 



dx 



0. (3) 



x=xb 

The so-called radiation boundary condition [SI E] is given by 



dp{x, t) 



dx 



= ^PixB,t). (4) 

x=xb 



This can be interpreted as a boundary which displays the characteristics of 
both the totally absorbing and refiecting cases i.e. there is partial absorption. 

In principle the diffusion equation can be solved with the radiation bound- 
ary condition by using the appropriate Green function, but it may be that 
this Green function is difficult to find. It then may be more convenient to use 
what is known as a ^sink term' in the Master equation instead. In general 
terms, a sink term represents the loss of probability density from a region of 
space as a consequence of absorption, attenuation, or some other similarly 
natured process. 

To see the equivalence explicitly jTU] we consider an absorbing region of 
width 2xb centered at the origin of the real line, at the boundaries of which 
an incident diffusive particle is absorbed with rate a. The Master equation 
for the probability density of the particle in this system reads 

dp{x,t) d^p{x,t) 

— — = — adixB - \x\)p{x,t) . (5) 

We define the survival probability of the particle at time t as the integral 
over all space 

r+oo 

n{t) = I p{x, t)dx (6) 



and probability that the particle is at a position |x| > -|- e as 

p{x,t)dx+ / p{x,t)dx . (7) 

-oo JxB+t 
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We assume that far away from the absorbing region that 



dp{x, t) 
dx 
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(10) 



(9) 



from ([T]). Now we demand that 




dn 



(11) 



to find 



dp 
dx 



x=±xb 



±—p{±XB,t) , 



(12) 



which is as the radiation boundary condition quoted above. Thus we have 
shown that a sink term of the form 



in a Master equation of the form of (fTl) is equivalent to the boundary condition 



As is comprehensively discussed in [7|, the length D/a is equivalent to 
the attenuation length in a composite medium. This is a medium in which 
there is a boundary at x = 0, and in the region x > particles can undergo 
free diffusion, and for x < the density of particles is attenuated [8]. It can 
be shown for this system that the probability density within the attenuating 
region decreases linearly to at the depth x = —D/a = —I, hence we refer 
to I as the attenuation depth. This furnishes a physical interpretation of the 
characteristic length scale D/a. 

2.2 Dimensionless variables 

We first define the inverse length scale 



aS{x — xb)p{x, t) 



(13) 




(14) 
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where ag ^ represents the characteristic diffusion length between resetting 
events. 

For convenience we summarise here some dimensionless combinations of 
parameters which will appear later: 



9 = 1 / — |xo| = «o|a;o| (15) 




(16) 

The dimensionless quantity 6 measures the ratio of the distance of the 
resetting position xq from the target (in this case at the origin) to the in- 
verse length ao- The dimensionless quantity (po measures the ratio of the 
attenuation length D/a to the characteristic diffusion length I/oq. 



3 Survival Probability Calculation 

To solve the Master equation ^ for the survival probability with partial 
absorption we perform the Laplace transform on the variable t. We define 
the Laplace Transform of the survival probability as 

POO 

q{z,s)= / g(z,t)e-^*dt . (17) 

We use this definition and the initial condition q{z, 0) = 1 to take the Laplace 
Transform of equation ([2]): 

We first find the solution to the homogeneous equation 

^^^0^-(r + .)g(^,.) = O, (19) 

which yields 

q{z,s) = Ae^' + Be-"' (20) 



D ' - (r + s)q{z, s) = -I - rq{xo, s) + aq{0, s)6{z) . (18) 



where 



«(^)=(^V^'. (21) 



(In particular from (14) we have ao = a{0).) From the full equation (18) (for 
2; 7^ 0) we find the particular solution 

q(z, s) = Ae"^ + Ee-"^ + ^^""^^^"'^^ . (22) 

r + s 
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We require that our solution be finite in the hmits z — ?■ ±00, which imphes 
separate solutions q{z, s)+ and q{z, s)_ for the regions where z is greater than 
and less than zero respectively. Using the required continuity at z = to 
find that A = 5, we obtain the result 



1 + rq{xo, s) 
r + s 



(23) 



Now we can use the discontinuity in the derivative which comes from the 
delta function in equation ^ to find A: 



D 



dq+ dq- 
dz dz 



ag(0, s) 



(24) 



2=0 



From (23) and (24) we derive the solution 

l-Q{z,s) 



where 



(s) is given by 



q{z,s) -- 
Q{z,s) 



s + rQ{xo, s) 

e-"kl 
" 1 + 2(j){s) ' 



(25) 
(26) 

(27) 



and 0(0) 



as given in (16). The long time behaviour of the survival 



probability can be determined from equation (25). As illustrated in Fig. [T| 
the analytic structure of this function in the complex s plane is a branch 
point s = —r, and a simple pole at s = sq which satisfies: 



where 



so(20(so) + l)+re-^^l^l = 



- r < So < . 



(2^ 



(29) 



In the large time limit, the residue from sq dominates the solution to the 
Bromwich inversion formula for this function. For convenience, we make a 
change of variable 

So = r{uo - 1) . (30) 



Under this substitution equation ( 28 ) becomes 



{uq - 1)(20oa/wo + 1) + exp 



D ' 



(31) 
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c 



branch cut 



branch point 



■ simple pole 



Re(s) 



Figure 1: Sketch of the singularities of the function (25) in the complex s 
plane. 



where 0o is given in (16). For further convenience we consider the case where 



z = xq (i.e. the initial position coincides with the resetting position) to find 



q{xo,t) 



(3mo - 1)00 + O(t>ouj\uo - 1) + 2n;/^ + e{uo - 1) 



1/2 



(32) 



The important point to note is that the survival probability decays expo- 
nentially and the rate of decay increases with absorption velocity a. If one 
expands to leading order in 1/a one finds 



Mo ~ M + 



4m*(1 



2{u*y/'^ - 9{1 - u*) a ' 



(33) 



where u* is the solution of (31) as a — )■ oo. 



4 Mean Time to Absorption 

The Mean Time to Absorption (MTA) T{z) for a particle which originated 
at 2; is [7] 



Tiz) 



dt 



q{z, t)dt = q{z, s = 0). 



(34) 



In our analysis we are considering mean first time to absorption of the 
searcher by the target, not just the mean first time to coincidence as in the 
case where the absorption is perfect. This means that it is possible for the 
searcher to pass through the target site without interacting with it. The 



MTA as defined above in (34) is still the appropriate measure of this process. 
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Using equation (25) and the definition in equation (34) we can write down 
an expression for tlie MTA: 



and in particular 



pQ!o|a;o| 

T{z) = [2(t>, + 1 - e-"ol^l] , (35) 



1 , , Opao\xo\ 

T{xo) = - [e""l^«l - 1] + 00 . (36) 



Tlius tlie effect of partial absorption is to increase the mean time to absorp- 



tion through the second term in (36). To understand better this term it 



is instructive to consider a more general problem which involves a resetting 
distribution^ V{x), i.e. the particle resets with rate r to a random posi- 
tion drawn from the distribution V{x). The resetting term rq{xo,t) in the 
backward Master equation ^ becomes r J dxV{x)q{x,t) so, for the survival 
probability q{z,t), the Master equation itself reads 

^^W^ = D^^^ - rq{z, t)+r j V{x)q{x, t)dx - ag(0, t)6{z) . (37) 

The calculation of the mean time to absorption is a straightforward general- 
isation of that presented above and the result is 



where 

y j dx'P(a;')e-"«l"-"'l 
p*{x) is the stationary distribution of the diffusive process with resetting 



p*[x) = y / dx'P(a;')e-"«l"-" I . (39) 



but without any absorption [3]. Thus (38) tells us that the effect of partial 



absorption is to increase the mean time to absorption by an additional term 
proportional to 1/a. In fact the precise form of the additional term is simply 
l/ap*(0), where p*(0) is the stationary probability of the searcher being at 
the origin (target site) in the presence of resetting but absence of absorption. 

4.1 Minimisation of the MTA with respect to resetting 
rate r 



Returning to (36 ) we are interested in finding the choice of r which minimises 



T[z). We start by considering the minimum of T{xo) with respect to r. 



9 



Setting = yields 

or 



(40) 



where 6 and 0o are defined in (15) and (16). The solutions of (40) are shown 
in Fig. |2j 



In the following we analyse the solution of (40) in the regimes a ^ a^D 



and a <^ a^D corresponding to strong and weak absorption. 



CD 




Figure 2: Plot of the value of 6 that minimises the MTA for a given (pQ. The 
value of 9 dX (po = 0, found numerically to be 9* = 1.5936, is indicated. 

First we note that for 0o ^ 1 (i.e. strong absorption, a ^ aoD) Eq. 



(40) reduces to the transcendental equation found in the study of total 
absorption ^4j: 

n* 

y = 1 - e-^* . (41) 
We then find that, to second order in 0o, 

of)* 

9 = 9*-2<Po + (42) 
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Figure 3: The minimised MTA, T*, as a function of (po, in units of x^/D. T* 



is approximately linear for both small and large 0o- (See (44) and (47). 



satisfies (40). It is helpful to rewrite Eq. (36) as 



T{xo) 



^0 



[(200 + l)e' - 1] , 



(43) 



and then, using Eq. (42 ), it can be shown that for 0o ^ 1 the MTA minimised 
with respect to r is 



T*(xo) 



Xn 



1 



■ 4 ^ 4(3 - e*) . 
1 + ^*^0 + 



(44) 



D e*{2-e*] 

To study the regime 0o ^ 1 (i-e- weak absorption, a <^ cto-D), we first rewrite 



Eq. (40) as 



1-9 



1 



1 



(45) 



^0 V 2^ 

Then we can see that for 0o 3> 1, — t- 1^, and explicitly to first order in 



11 



1/00 

1 — 2e~^ 

^ ^ 1 + — — . (46) 
200 

Using this result we find that the minimised MTA is, to leading order in 

1/00, 

T*(xo)^2e^0o. (47) 
We plot T* as a function of 0o in Fig. [3j 



5 Many Independent Searchers 

We now consider the multiparticle version of the search process comprising 
a single immobile target at the origin and many searchers which are initially 
uniformly distributed on the line with uniform density p. The searchers are 
independent of each other and the position of each searcher evolves stochas- 
tically starting at its own initial position to which it resets. 
The survival probability of the target, Ps(t), is given |4j by 

TV 

Ps{t) = l[q{x„t) , (48) 
1=1 

where q{xi,t) is the survival probability in the single searcher problem. The 
positions Xi are independent and distributed uniformly within the box [-L/2, 
L/2]. The average probability Pg^it) = {Ps{t))x and the typical probability 
Pl^^it) = exp[{\rLPs(t))x], where {■)x denotes averages over Xj's. 

In random additive processes the average of the random variable and its 
typical or most probable value exhibit the same behaviour. For multiplica- 
tive processes (where a product of random variables is considered such as 
is described here), there exist extreme events which, although exponentially 
rare, are exponentially different from the typical value of the product [12]. 
Thus, for multiplicative processes, it is important to consider both average 
and typical values. 

5.1 Average Survival Probability of the Target 

The average survival probability of the target, P^'"{t), is given by 

Pnt) = (g(x,t))f = exp(iVln[l - (1 - g).]) , (49) 

where 

(l-^)x = 7/ dx[l-qix,t)]. (50) 

J-L/2 
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Taking N, L oo while keeping the density of searchers p = N/ L fixed 
(using q{x,t) = q{—x,t)) we find: 



P^%t) exp ^-2p ^ dx[l - q{x, t)]^ = exp( 



-2pM(t)) . (51) 



Defining 

POO 

M{s) = / M{t)e-''dt , (52) 



we can use equation (25) to find 

M(s) = \n(l + . (53) 

^ ' sra V s(20(s) + 1)) ^ ' 

In principle this may be inverted to find P^'"{t). Rather than present the 
formula which takes the form of a double convolution it is more instructive 
to examine the effect of partial absorption on asymptotic behaviour directly 
from (53). 

For convenience we first write 

0(s) = k{r + s)2, where k = , (54) 



so that we can then expand (53) to leading order to find 

M(s) = W — - (-Ins + In r -ln[2A;v^+ 1] +0(s)) , (55) 
V r s 



and using the identity [13] 

^_,r Ins 



lnt + 7, (56) 

s 

where 7 is Euler's constant, we find then that for long times 

P"'^(t) ~ exp [-2p (Inrt + 7 - In (20o + 1))] • (57) 

Thus the effect of partial absorption is to change the asymptotic decay of the 
average probability by a multiplicative factor (1 + 20o)^''- 
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5.2 Typical Survival Probability 

The typical survival probability of the target, P^y^it), can be expressed as 



N 



Pi^^(t)=exp5^(lng(x, , t))^ = exp 



i=l 



rL/2 

2p / dx\n q{x,t) 
Jo 



(5^ 



In the long-time limit we can use Eq. (32) to show that 



/ dxo In g(xo, t) — const — t / dxo|so(xo)| 
Jo Jo 



(59) 



and P^yp{t) ~exp(-2/pt). 

We can calculate the integral / = dxoSo{xo) using the following pro- 
cedure. First, using equation (28) we make the substitution z = —so/r to 
find that z satisfies 



z (0o(l - zY^' + 1) = exp(-0(l 



which we use to rewrite J, by eliminating 6, as an integral over dz: 



where 





^rD I dz 



Inz ^ln[(0o(l-^)^/' + l)] 



(1 - z)l/2 
1 



(1 - z)l/2 



We make the substitution 
to obtain 



y=i-z 



(60) 

(61) 

(62) 
(63) 



/yo 
dy (ln(l - 2/2) + in[(0o2/ + 1)]) , where = 1 - Zq . (64) 

The first term in this integral, which will henceforth be labelled /q, 

/yo 
dy [ln(l -y) + ln(l + y)] (65) 



which can be evaluated to obtain 

;i + yo) 



4VrD 



1 -ln2-?/o + 



ln(l + yo)-^^^^ln(l-2/o) 



(66) 
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For 00 ^ 1 we find 



(67) 



and consequently 



(6^ 



To calculate the second term, Ji, in (64) we again use the substitution 



(|63J) to find 



h = [(00 + 1) In (00 + 1) - (2/000 + 1) In (0ol/o + 1) - 0o (1 - 2/o)] • 



(69) 



From (67) we see that, for 0o ^ 1, 



Ji ~ -AVrD 



(70) 



Combining the results (68) and (70) we find that for 0o ^ 1 



I ~ AVrD 1 - In 2 - 



(71) 



In this limit, from (68) and (70), the correction to the total absorption decay 
rate is linear in 00 : 



P'yP{t) ~ exp 
The full result for P^yp{t) is 



-8VrD ( 1 - ln2 - ) pt 



(72) 



p'yp{t) 



exp — 8ptv rD 
6o + l 



l_ln2-yo + ^^^ln(l + l/o) 



:i - yo) 



In (00 + i: 



yO0O + 1 



In (0o?/o + 1) - 



2 

1 - yo 



ln(l 



6 Conclusion 



In this work we have considered the dynamics of a diffusive searcher in a 
system with a partially absorbing target, as defined in ([TJ, which gives a 
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more realistic description of many, varied search processes. Our study has 
revealed some straightforward, significant and intuitive consequences for the 
dynamics as a direct result of imperfection in the absorption at the target. 

We see from Section |3| and equations (32) and (33) in particular, that 
the survival probability of the searcher (or target) decreases exponentially 
with time, with a decay constant which increases as the absorption constant 
a increases. As the target comes closer to being perfectly absorbing, the 
survival probability decays away much sooner. 

Our study of the MTA in Section |4] has revealed that the mean time to 
absorption is increased by an additive term l/ap*(0) (see (38)). As a de- 
creases and the target is made more 'imperfect' the mean time to absorption 
increases, as expected. 

For multiplicative processes, such as the many searcher problem analysed 
in Section [5} the distinction between typical and average probabilities be- 
comes significant. This is emphasised by the difference in the form of effects 
of imperfect absorption on the typical and average survival probabilities of 
the target in the many searcher system. We see from Section 5.1 that P"''"{t) 
is modified by multiplicative factor (1 + 20o)^''; whereas in Section 5.2 it is 



the decay rate of P*^^(t) which is decreased by factor proportional to 1/a. 

An important quantity that emerges from this this work is the dimension- 
less ratio 00 = cto-D / a. This quantity is a ratio of length scales characteristic 
to the system: the length scale 1/ao = \/ D jr has already been established 
as the characteristic displacement of the searcher due to diffusion between 
reset events|3j; the length D/a is the attenuation depth in the composite 
medium discussed in [7]. Thus the dimensionless variable 0o represents the 
ratio of these two length scales, and the dynamics are only modified in terms 
of this ratio and not the absolute strength of the imperfection which has 
been introduced to the target. It is the competition between resetting and 
the absorption that controls the dynamics of the system. 

In future work it would be of interest to extend the present calculations 
to higher dimensions. Also a related problem of interest [2] would be to con- 
sider discrete time jump processes where the searching particle may jump 
over the target thus reducing the effective absorption rate. 
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